Thickness-Magnetic Field Phase Diagram at the Superconductor-Insulator Transition 
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The superconductor-insulator transition in ultrathin films 
of amorphous Bi was tuned by changing the film thickness, 
with and without an applied magnetic field. The first experi- 
mentally obtained phase diagram is mapped as a function of 
thickness and magnetic field in the T=0 limit. A finite size 
scaling analysis has been carried out to determine the critical 
exponent product i/z, which was found to be 1.2 ± 0.2 for the 
zero field transition, and 1.4 ± 0.2 for the finite field tran- 
sition. Both results are different from the exponents found 
for the magnetic field tuned transition in the same system, 
0.7 ±0.2. 

PACS numbers: 74.76.-w, 74.40.+k, 74.25.Dw, 72.15.Rn 



Superconductor-insulator (SI) transition in ultrathin 
films of metals is believed to be a continuous quantum 
phase transition [Q which can be traversed by changing 
a parameter such as disorder, film thickness, carrier con- 
centration or the applied magnetic field [§,[§. The scal- 
ing theory and a phase diagram for a two-dimensional 
system as a function of disorder and magnetic field was 
postulated by Fisher et al. , based on the assumption 
that this transition can be fully described in terms of a 
model of interacting bosons, moving in the presence of 
disorder. The dirty boson problem has been extensively 
studied by quantum Monte Carlo simulations p j^] , real- 
space renormalization group calculations pjJlC |, strong 

but 



coupling expansion Jll| and in other ways ]12|-|1 Z 
there is still some disagreement as to the universality 
class of the transition. Conflicting experimental evidence 
suggests that the bosonic model might be relevant [ fL5| , 
but does not give the full picture An alternative 

model of interacting electrons has also been proposed 
p7[ . Experimentally, the thickness tuned transition has 
been studied in the context of the scaling theory in zero 
magnetic field |l^] . In the present work, for the first time, 
the SI transition was tuned by systematically changing 
the film thickness in a finite magnetic field. This allows 
us to map a phase diagram as a function of thickness 
and magnetic field in the T=0 limit and to determine 
the critical exponents using a finite size scaling analysis 
at different fields. The results suggest that this transition 
is similar to the zero field transition, but the exponent is 
different from that of the magnetic-field tuned transition 
studied on the same set of films |ll| . 

The ultrathin Bi films were evaporated on top of a 10 



A thick layer of amorphous Ge, which was pre-deposited 
onto a 0.75 mm thick single-crystal of SrTiO^ (100). The 
substrate temperature was kept below 20 K during all 
depositions and all the films were grown in situ under 
UHV conditions (~ 10~ 10 Torr). Under such circum- 
stances, successive depositions can be carried out with- 
out contamination to increase the film thickness grad- 
ually in increments of ~ 0.2A Film thicknesses were 
determined using a previously calibrated quartz crystal 
monitor. Films prepared in this manner are believed to 
be homogeneous, since it has been found that they be- 
come connected at an average thickness on the order of 
one monolayer |pOf . Resistance measurements were car- 
ried out between the depositions using a standard dc 
four-probe technique with currents up to 50 nA. Mag- 
netic fields up to 12 kG perpendicular to the plane of the 
sample were applied using a superconducting split-coil 
magnet. 

The evolution of the temperature dependence of the 
resistance as the film thickness changes is shown on Fig. 
[l]. The thinnest films show an exponential temperature 
dependence of the resistance at low temperatures, consis- 
tent with variable range hopping, which crosses over to 
a logarithmic behavior for thicker films p3j. At some 
critical thickness, d c , the resistance is independent of 
temperature, while for even thicker films it decreases 
rapidly with decreasing temperature, indicating the on- 
set of superconductivity. The critical thickness can be 
determined by plotting the resistance as a function of 
thickness for different temperatures (inset of Fig. |l|) and 
identifying the crossing point for which the resistance 
is temperature independent, or by plotting dR/dT as a 
function of thickness at the lowest temperatures and find- 
ing the thickness for which (dR/dT) = 0. 

In the zero temperature quantum critical regime the 
resistance of a two dimensional system is expected to 
obey the following scaling law [IJQ]: 

R(S,T) = R C fiST- 1 /^) (1) 

Here 6 = d — d c is the deviation from the critical thick- 
ness, R c is the critical resistance at d = d Cl f(x) is a 
universal scaling function such that /(0) = 1, v is the co- 
herence length exponent, and z is the dynamical critical 
exponent. We rewrite Eq. |l| as R(S,t) = R c f(St), where 
t = r £^ 1 l vz , and treat the parameter t(T) as an unknown 
variable which is determined at each temperature to ob- 
tain the best collapse of all the data. The exponent vz is 
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then found from the temperature dependence of t, which 
must be a power law in temperature for the procedure 
to make sense. This scaling procedure does not require 
detailed knowledge of the functional form of the temper- 
ature or thickness dependence of the resistance, or prior 
knowledge of the critical exponents. It is simply based on 
the data which includes an independent determination of 
d c . 

The collapse of the resistance data as a function of 
St in zero field is shown in Fig. ^. The critical expo- 
nent product vz, determined from the temperature de- 
pendence of the parameter t (inset of Fig. |^), is found 
to be vz = 1.2 ± 0.2. This result is in agreement with 
the predictions of Ref. from which z=l would be 

expected for a bosonic system with long range Coulomb 
interactions independent of the dimensionality, and v > 1 
in two dimensions for any transition which can be tuned 
by changing the strength of the disorder A simi- 

lar scaling behavior has been found in ultrathin films of 
Bi by Liu et al. []l8|], with the critical exponent product 
vz Ki 2.8 on the insulating side and vz ps 1.4 on the su- 
perconducting side of the transition. The fact that vz 
was found to be different on the two sides of the transi- 
tion raises the question of whether the experiment really 
probed the quantum critical regime. We believe that the 
scaling was carried out too deep into the insulating side, 
forcing the scaling form (Eq. [I]) on films which were in a 
fundamentally different insulating regime pl[ . Such films 
should not be expected to scale together with the super- 
conducting films, hence the discrepancy on the insulating 
side of the transition. In the present work, the mea- 
surements were carried out at lower temperatures than 
previously studied and with more detail in the range of 
thicknesses close to the transition. We were able to scale 
both sides of the transition with vz ss 1.2, which is close 
to the value obtained by Liu et al. on the superconduct- 
ing side of the transition. Our result is also in very good 
agreement with the renormalization group calculations 
P, [iol|b!f| , and close to that found in Monte Carlo sim- 
ulations by Cha and Girvin Jjj, S0rensen et al. j(| and 
Makivic et al. J§. 

In addition to the above, the magnetoresistance as a 
function of temperature and magnetic field was measured 
for each film. By sorting the magnetoresistance data, one 
can probe the thickness-tuned superconductor-insulator 
transition in a finite magnetic field, which has not been 
studied before. The same analysis as described above 
was carried out for several fixed magnetic fields, rang- 
ing from 0.5 kG up to 7 kG. The normalized resistance 
data as a function of the scaling variable for six different 
values of the magnetic field shown on Fig.^ all collapse 
on a single curve, which suggests that the scaling func- 
tion is indeed universal. The critical exponent product 
determined from the temperature dependence of the pa- 
rameter t (inset of Fig. ||) is found to be vz = 1.4 ± 0.2, 
apparently independent of the magnetic field. An applied 



magnetic field is generally expected to change the uni- 
versality class of the transition, since it breaks the time 
reversal symmetry. We find, however, that the critical ex- 
ponent product vz for the thickness driven SI transition 
in a finite magnetic field is very close to that obtained 
for a zero-field transition, given the experimental uncer- 
tainties. This result is in agreement with Monte Carlo 
simulations of the (2+l)-dimensional classical XY model 
with disorder by Cha and Girvin [Q, which find vz « 1.07 
for the zero-field transition and vz 1.14 in a finite mag- 
netic field. The critical resistance R c at the transition is 
non-universal, as it decreases with increasing magnetic 
field. 

Furthermore, once a magnetic field is applied and 
the time-reversal symmetry broken, the thickness-tuned 
transition in the field is expected to be in the same uni- 
versality class as the transition which is tuned by chang- 
ing the magnetic field at a fixed film thickness. The 
magnetic-field tuned transition was studied on the same 
set of films used in this study which allows a di- 
rect comparison of the critical exponents, without hav- 
ing to worry about differences in the microstructure be- 
tween different samples. The resistance as a function 
of temperature for five selected films from Fig. [I] was 
studied in magnetic fields up to 12kG applied perpen- 
dicular to the plane of the sample. The critical expo- 
nent product was determined using the method described 
above, but with the magnetic field as the tuning parame- 
ter rather than the film thickness. It was found to be 
vb z b = 0.7 ± 0.2, independent of the film thickness, 
which strongly suggests a universality class different from 
that of the thickness-tuned transition, both in zero-field 
and in the field. This result does not agree with the pre- 
dictions based on the model of interacting bosons in the 
presence of disorder (30. It does, however, agree with 
what might be expected from a similar model without 
disorder ||,[7j . The details of the analysis and this unex- 
pected value of the critical exponent product, as well as 
the discussion of its disagreement with previous determi- 
nations [^3|, are reported elsewhere |l9 |. 

The fact that the SI transition was traversed by chang- 
ing both the thickness and the magnetic field indepen- 
dently on the same set of films allows us to determine 
the phase diagram as a function of thickness and the 
magnetic field in T=0 limit, which is shown in Fig. ^. 
The films characterized by parameters which lie above 
the phase boundary are "superconducting" (SR/5T < 
at finite temperatures), and the ones bellow it are "insu- 
lating" (SR/ST > at finite temperatures). The phase 
boundary itself follows a power law: d c (B) —d c (0) oc B*, 
where x ~ 1.4. Our results imply that the critical expo- 
nent product vz depends on whether this phase bound- 
ary is crossed vertically (changing the thickness at a con- 
stant magnetic field), in which case vz ~ 1.4, or horizon- 
tally (changing the magnetic field at a fixed thickness), 
in which case vbzb ~ 0.7. In other words, relatively 
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weak magnetic fields which are experimentally accessible 
to us do not significantly change the universality class 
of the thickness-tuned transition, but the magnetic-field 
tuned transition is in a different universality class from 
the thickness-tuned transition in B ^ 0. In the absence 
of a detailed theory, we can only speculate on the origins 
of this surprising result. 

The first important issue we wish to discuss is the 
role of the film thickness as the control parameter. 
Adding metal sequentially to a quench-condenscd film 
has been shown to decrease the disorder, since the in- 
creased screening smooths the random potential seen by 
the electrons. It presumably also increases the carrier 
concentration, which in the presence of an attractive 
electron-electron interaction might result in an increased 
Cooper pair density. 

Increasing the film thickness might therefore be 
thought of as adding Cooper pairs, which condense at 
some critical density. In a similar way, an applied mag- 
netic field adds vortices, which behave as point parti- 
cles and also condense at some critical density, making 
the system insulating. However, this symmetry between 
charges and vortices is not perfect, since Cooper pairs 
interact as 1 /r and vortices interact logarithmically pl| . 

If the mechanism responsible for the localization in the 
magnetic-field tuned transition is different from that of 
the thickness-tuned transition, than having a non-zero 
magnetic field may not play a major role in the thickness- 
tuned transition. Also, the correlation length associated 
with the thickness-tuned transition would then be differ- 
ent from that associated with the magnetic-field tuned 
transition. The disorder might be important in one case 
and not in the other, depending on how these correlation 
lengths compare to the length scale which characterizes 
the disorder. 

A second difference between the thickness-tuned and 
the field-tuned transitions may be the nature of disor- 
der itself. In the field-tuned transition the geometry of 
the film is fixed, and the disorder does not change. In 
the disorder-tuned transition, each film in the sequence 
of films will have a slightly different microstructure, so 
that the disorder may have to be averaged over the dif- 
ferent configurations. It has been suggested recently ]25j] 
that the nature of the disorder averaging might play an 
important role in determining the critical exponents. 

Another possibility is that the localization transition is 
enhanced by percolation effects |2q] as the film thickness 
is tuned. This approach takes into account local fluc- 
tuations of the amplitude of the superconducting order 
parameter, which are routinely neglected in the scaling 
theory and the numerical simulations. Percolation of is- 
lands with strong amplitude fluctuations might change 
the localization exponents obtained from the scaling the- 
ory |^6| . The role of percolation has also been empha- 
sized in recent discussions of low temperature transport 
in these systems p7| . 



Finally, one cannot exclude the possibility that to ac- 
cess the quantum critical region these measurements need 
to be carried out at much lower temperatures or at high 
frequencies p8| |. 
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FIG. 1. Resistance per square as a function of temperature 
for a series of bismuth films with thicknesses ranging from 
9A (top) to 15A (bottom). Inset: Resistance as a function of 
thickness for the same set of films close to the transition at 
low temperature. Different curves represent different temper- 
ature, ranging from 0.14K to 0.40K. 

FIG. 2. Resistance per square as a function of the scaling 
variable, t|d-d c |, for seventeen different temperature, ranging 
from 0.14K to 0.5K. Here t = T~ 1/vz is treated as an ad- 
justable parameter to obtain the best collapse of the data. 
Different symbols represent different temperature. Inset: Fit- 
ting the temperature dependence of the parameter t to a 
power law determines the value of vz. 

FIG. 3. Normalized resistance per square as a function of 
the scaling variable, t\d — d c \, at different temperature, rang- 
ing from 0.14K to 0.5K. Different symbols represent different 
magnetic fields, ranging from 0.5kG-7kG. Here t = T~ x l" z j s 
treated as an adjustable parameter to obtain the best collapse 
of the data, and R c is the resistance at d = d c . Inset: Fitting 
the temperature dependence of the parameter t to a power 
law determines the value of vz. 



FIG. 4. The phase diagram in the d-B plane in the T=0 
limit. The points on the phase boundary were obtained 
from disorder driven transitions (triangles) and magnetic field 
driven transitions (circles). The solid line is a power law fit. 
The values of the critical exponent product are shown next 
to the arrows giving the direction in which the boundary was 
crossed. Here d c is taken to be the critical thickness in zero 
field. 
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